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Most of us will have at some time thrown a pebble into water, and watched the ripples spread 
outwards and fade away. But now there is a way to reverse the process, and make those ripples 
turn around and reconverge again, ... and again, and again. To do this we have designed the 
Maxwell's Fishpond, a water wave or "Transformation Aquatics" version of the Maxwell's Fisheye 
lens [[j], ^] that is now well-known from transformation optics. These are transformation devices 
where wave propagation on the surface of a sphere is not modelled on an actual sphere, but in 
a flat device with carefully designed spatially varying properties. And just as for rays from a 
point source on a sphere, a wave disturbance in a Maxwell's Fisheye or Fishpond spreads out at 
first, but then reforms itself at its opposing point in the device. Here we show how such a device 
can be made for water waves, partly in friendly competition with comparable electromagnetic 
devices [[| and partly as an accessible and fun demonstration of the power of transformation 
mechanics. To the eye, our Maxwell's Fishpond was capable of reforming a disturbance up to 
five times, although such a feat required taking considerable care, close observation, and a little luck. 



What can you see in the video at [ittp:/ /www. qols.ph.ic.ac.uk/~kinsle/files/MFishpond/ 



I. INTRODUCTION 

Transformation Optics and Transformation Acoustics 
are powerful new techniques used to design transforma- 
tion devices, which usually tend to the exotic - invisibility 
cloaks, illusion generators, and so on. Regardless of their 
utility, none seem to be useful as demonstration devices 
accessible to non-specialists. 

Here, in contrast to this typical situation, we describe 
how to make a transformation device (T-device) that fits 
on a tabletop and which controls water waves visible to 
the naked eye. Although most T-devices such as cloaks 
cannot be made with simple isotropic materials, one im- 
portant type can: one that transforms from waves trav- 
elling on the surface of a uniform sphere, by squashing 
them down onto waves travelling on a disk. This trans- 
formative squashing can preserve the properties of the 
original wave propagation - that of circular trajectories 
of equal circumference - by making the disk properties 
non- uniform: here, we make a shallow pond of appropri- 
ately varying depth. 

In optics, such a device is known as the Maxwell's Fish- 
eye lens |2|] , which has a long history, originally being 
proposed in a problem set in the 1850's Q |). Until re- 
cently only an obscure theoretical curiosity, the device 
was brought to wider attention by controversial claims 
that in EM such a device could generate perfect images 
[||. Leaving aside that debate (see |7j for a recent crit- 
ical summary), Leonhardt usefully noted that actually 
building the device becomes possible if you take only the 
central portion and surround it by a mirror [|| . This trick 
also makes the water wave version easier to build. But 
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what makes a Maxwell's Fisheye lens - or its Fishpond 
counterpart - interesting? 

As in an ordinary pond, a pointlike wave source any- 
where on the surface of a Maxwell's Fishpond generates 
an outgoing set of ripples. However, in a Maxwell's Fish- 
pond, irrespective of where they were generated, the rip- 
ples do not just spread out and disperse, they instead re- 
converge on the opposite side, before again diverging and 
travelling back to the start where they again reconverge, 
and repeating this process until they eventually dissipate 
- just as they do in the optical Maxwell's Fisheye lens 
This matches exactly what rays or waves confined on the 
surface of a sphere would do, which as described above 
is the design principle behind both devices. More gener- 
ally, the Maxwell's Fisheye lens is one of a moregeneral 
class of classical transformation optics devices U; and, 
as we have done here, these others, such as the Eaton or 
Luneberg lens, will also have accessible transformation 
aquatics equivalents based on water waves - just as they 
do on the more technologically exotic field of plasmonics 

§0- 

Here we will derive the water depth profile needed to 
make a Maxwell's Fishpond in the shallow water limit, 
and compare that profile to the simple approximation 
of a shallow spherical dome. Although an approximate 
spherically bottomed pond can work surprisingly well, 
our accurate device can hint at - at least to the eye - 
up to five successive refocussings! More rigorously, we 
also present simulation results indicating how the device 
works in practise, as well as two experimental schemes 
set up with relatively litte demands on equipment. 

Contributions to this work were as follows: PK con- 
ceived the Maxwell's Fishpond idea, and built the first 
crude prototype. He also designed the version used here, 
but with students NK and TT shadowing that design 
process. NK and TT did the first experiments, and CT 



1 



FISHPOND 



Maxwell's Fishpond 



[Dr. Paul. Kinsler@physics.org" 
http://www.kinsler.org/physics/ 



and JT followed next; all four writing reports and giv- 
ing presentations as part of their coursework. PK was 
the primary author of this paper, assisted by material 
from the student reports, he also did all the computer 
simulation work. CT, JT, and TT also assisted in the 
preparation of this final manuscript. 



II. FISHEYE, FISHPOND 

The Maxwell's Fisheye concept is based on mimick- 
ing the properties of ray trajectories on the surface of 
a sphere using a flat surface with spatially modulated 
properties. This is interesting, because on a sphere all 
geodesies are the so-called "great circles" , and any set of 
rays emitted from a point follow these and then automat- 
ically re-focus on the exact opposite side of the sphere. 
Thus, any flat version would also have this property - 
rays leaving from any point would automatically refo- 
cus at the complementary point of the plane. Thus, in 
the fisheye, an object at any point is guaranteed to form 
an image; this is most certainly not the case in ordi- 
nary imaging systems. Further, the rays would then re- 
diverge, only to reconverge again; in an ideal ray device, 
this would continue forever. 

To achieve the transformation from a spherical device 
to a flat one, we use a stereographic projection. Imagine 
the sphere sitting with the south pole on a flat sheet, as 
shown on the upper part of fig. [l|. Then each point A or 
B on the sphere is mapped onto the sheet by following 
a straight line from the north pole, through A (or B), 
and onward until it intersects the sheet at A' (or B'). 
In this way the southern hemisphere maps onto a disk 
on the sheet centered on the south pole; the northern 
hemisphere is mapped to points further away; with points 
very near the north pole being extremely remote, and the 
north pole itself having to be omitted. 

Of course, although we would like to make a fisheye (or 
fishpond) based on this projection, we do not want one 
that is infinitely big. To solve this conundrum we follow 
Leonhardt Q and place a mirror at the equator, confining 
all ray paths to the southern hemisphere, and so confin- 
ing all projected rays to inside a circle with twice the 
sphere's radius. Since both hemispheres have the same 
properties, the ray properties are preserved - although 
the great circles are now folded back on themselves and 
have a kink where they are reflected, they still are guar- 
anteed to form a image of any point. 

The process of opening out and flattening the surface 
of a sphere into an equivalent sheet, as if it were a map 
projection used for an atlas, has an important feature. 
Regions near the equator are stretched and expanded, 
while those near the south pole are only slightly changed. 
Mathematically we can define a complex quantity z = 
x + ty, where x, y represent the Cartesian coordinate on 
the plane. This means that any point (X, Y, Z) (or at 
angles <p) on the unit sphere is projected (or "mapped") 



SIDE VIEW 
OF SPHERE 




9J9t[dSTUI9t[ msmjou 
jo iioipafoid 



FIG. 1: The sphere-to-plane fisheye projection. The south- 
ern hemisphere of the sphere of radius s becomes a finite disk 
of radius ro = 2s, and the northern hemisphere becomes the 
entire plane that remains outside the disk. Lines from pole- 
to-pole (meridians, or lines of longitude) on the sphere be- 
come radial lines (see e.g. the dot-dashed line), circles on the 
sphere parallel to the equator (parallels, or circles of latitude) 
become projected circles whose size depends on how far north 
or south of the equator (dashed line) they are. 



down onto 



X + iY 
l-Z 



= exp (i</>) cot (0/2) . (1) 



This means that a given line element dS 2 = dX 2 + dY 2 + 
dZ 2 on the sphere is transformed into a line element in 
the plane dR 2 = dx 2 + dy 2 that progressively lengthens 
as we move towards the equatorial perimeter at r = 1. 



1 



rdR 2 = dS 2 . 



(2) 



where r 2 = x 2 + y 2 , and we can also note that angles on 
the sphere are preserved when projected onto the plane. 

This length transformation means that an object (or 
ray) travelling at a fixed speed on the sphere, if projected 
onto the disk on the plane, would have to travel faster 
on the plane, the nearer it got to (the projection of) 
the equator. Thus a fixed object speed vq on the sphere 
is projected onto the disk as a radially varying velocity 
profile v(r). Because of the way the projection works 
- or the line element dS converts to dR - the velocity 
profile v(r) for a sphere of radius ro/2 is projected onto 
its counterpart disk of radius ro as 



v(r) v 



(3) 
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This velocity profile has a counterpart in a gradient- index 
version of Snell's law that steers propagating rays so that 
they match the paths that follow from projections of the 
great circle paths on the sphere. Indeed, from a mathe- 
matical perspective we might have expected this, because 
the projection preserves angles, which also means that 
the device can be made with isotropic materials. 



(a) 



200 mm 








7.5 mm 



A. Optics 

In optics, one can actually make a thin spherical shell 
that will guide light inside it (see e.g. [fTl|), just like the 
spherical reference device for a planar fisheye lens. To 
obtain a design for the Maxwell's Fisheye lens we need 
only work out how to design an optical device that has a 
light speed profile defined in eqn. (Jfy by modulating the 
refractive index. If starting with a shell with refractive 
index no (and hence speed of light d = c/tiq), we can 
convert this to a disk with a radially varying refractive 
index profile 



n(r) = 



no 



1 + (r/roY 



(4) 



where no is the maximum refractive index we can achieve, 
and ro is our desired radius scale. In the original fisheye 
lens, r was unbounded, causing the device to need un- 
realistically small values of n at large radii r, and even 
before then the inconvenient requirement for n < 1 at 
r > vq. The introduction of a mirror, as discussed above, 
circumvents this restriction; and for no > 2 the minimum 
refractive index required is always > 1. 

Electromagnetic Maxwell's Fisheye lenses have been 
made (e.g. || [l2]-[l4|), but require significant techno- 
logical skill to build and investigate. Hence our interest 
in water waves, which allow access to these interesting 
devices to a much wider audience. 



B. Water waves 

We want to make a Fishpond, not a Fisheye; water 
waves are easily visible, intuitive, low-tech, and are ac- 
cessible and safe for a wide variety of ordinary people. 
Nevertheless, experimental water wave systems can still 
be used to model - at least in part - quite a surprising 
variety of phenomena: e.g. event horizons and Hawking 
radiation 



16 and neutron star collapse 



To obtain a design for such a Maxwell's Fishpond we 
need only work out how to design a device that has the 
speed profile defined in eqn. (||). In general, water waves 
can have a complicated and nonlinear behaviour, so con- 
structing a general fishpond will be either very difficult 
or impossible. But there is an important subset of wa- 
ter waves for which we can get a simple solution - those 
waves that occur in very shallow water. 

For water of a constant depth that is significantly less 
than a wavelength, the wave speed for small waves is 




FIG. 2: The Maxwell's Fishpond, (a) A cross section with 
the vertical scale grossly exaggerated for clarity. The indi- 
cated dimensions are those used for our actual device, but 
heights or widths can be rescaled freely - subject to the pro- 
viso that ripples will have a wavelength longer than the maxi- 
mum depth. In fact, an even shallower fishpond would provide 
a better match to this criterion, but since water has a signif- 
icant surface tension, this makes covering the centre region 
problematic, (b) A photograph of our device. 



simply 



gd, 



(5) 



where g = 9.81m/s 2 is the gravitational acceleration and 
d is the water depth. In this extreme limit, other factors 
such as the wave amplitude and wavelength no longer 
matter, and we can control any suitable wave with the 
same depth modulation. Further, as long as d(r) varies 
slowly over wavelength scales, we can use the formula to 
describe waves travelling across a varying depth profile 
d(r). 

If the centre of the fishpond has depth do, the water 
wave speed there is ^(0) = \/gdo. Thus the radial wave 
velocity profile will be 



(6) 



v w (r) = v w (0)Vd, 



where d(r) = d(r)/do is the relative depth profile. 

Thus, comparing eqns. (||) and (||) we see that to 
match the two velocity profiles we need to have 



d(r) 
d(r) 



1 + (r/r y 
1 + (r/r y 



l + 2(r/r f + (r/r ) 4 



(7) 

(8) 
(9) 



The parameters we chose for our Fishpond were based on 
an assumed water wavelength of about 2cm; the result 
can be seen in fig. ||. 

One might also construct Eaton and Luneburg lenses Q 
(or even their generalizations Q ) for shallow water waves 
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based on a reference depth do and reference radius 7*0. By 
comparing velocity profiles between the optical and water 
wave cases, we find that the retro-reflecting Eaton pond 
feature and the focussing Luneburg pond feature need 
the depth profiles 

^EatonO) = ~ , (10) 

zro — r 

^Luneburg(r) = ^ ■ (11) 



C. Other waves 

It is possible to imagine other implementations of the 
Maxwell Fisheye concept. For example, an adaption of 
the expertise demonstrated by Bramhavar et al in |l8[ 
might give rise to an appropriately tapered "Maxwell's 
Platter" with the same behaviour for acoustic waves in a 
solid. 



III. MODELLING 

We tested our design using computer simulations for 
a variety of cases ranging from those applicable to the 
ideal Maxwell's Fishpond and approximate Fishponds, to 
a full finite element simulation for our Fishpond device. 
For the idealized comparisons, we used the fact that the 
fisheye and Fishpond behave in an essentially identical 
manner, once the distinctions between the polarizable 
EM field and scalar water waves have been accounted for. 
This means that since an ideal Fishpond has the same 
properties as an ideal fisheye lens, FDTD simulations of 
Maxwell's equations for the fisheye lens will indicate the 
behaviour of a suitably idealised Fishpond. 

To get an initial estimate of the importance of the cor- 
rect depth profile, we used MEEP |l§] FDTD simulations 
of Maxwell's equations. This approach was taken be- 
cause we already had such EM simulations running, and 
because the MEEP software is flexible, open source, and 
freely available. All that is required is to converted our 
chosen depth profiles back into a a refractive index pro- 
file using the reverse of the process that led to eqn. (^). 
As well as the exact Maxwell's Fishpond, and amongst 
other variations, we modelled an approximate depth pro- 
file based on a dome or section of a sphere. It turns out 
that as long as the correct 1:4 ratio of minimum to maxi- 
mum depths is maintained, this worked remarkably well. 

In fig. H we can see the simulation results equivalent 
to our shallow water wave model, showing snapshots that 
indicate the image reformation properties. We see that 
the ideal Maxwell's fisheye lens and Fishpond will give ac- 
curate refocussing (fig. |](a,b)), and this is repeated many 
times before the performance starts to degrade due to the 
dispersion caused by how different wavelengths interact 
with the finite-sized geometry, (see e.g. ||). Next, the 




FIG. 3: Snapshots from simulations representing an ideal 
Maxwell's Fishpond (a,b), an approximate Fishpond with a 
spherical (domed) depth profile (c,d), and an ordinary flat- 
bottomed fishpond (e,f). The upper frames (a, b) show the 
fisheye/pond wave patterns near the first and second refo- 
cussing times; the middle ones (c,d) the approximate SB 
pond, and the lower two (e,f) show the non-focussing be- 
haviour of a flat-bottomed pond when the ripples are not 
started at the exact centre. To the eye, the second Fishpond 
reformation (b) is essentially identical to the original source 
wave. 



simulations matching the approximate domed pond pro- 
file, do quite well (see fig. [|(c,d)), but with some distor- 
tion clearly evident on the second reformation. However, 
the chosen "best reformation" snapshots of the domed 
pond flatter slightly, as the frames before and after shown 
a significant ellipticity; and the third reformation (not 
shown), whilst still giving a localised wave bunch, has 
lost its concentric-ring character. Finally, simulations of 
a flat-bottomed pond (see fig. ||(e,f)) show only a poor 
attempt at a first focus, followed by a rapid evolution 
towards an apparently random pattern with no reforma- 
tions apparent at all. Other simulations including a vari- 
ety of strengths of non-radial distortion of the Fishpond 
depth profile were also performed. Depth variations of 
about 10% away from the exact profile do give tolera- 
ble results for the first few image reformations, but the 
distortion strongly degrades the beautiful concentric-ring 
character of the exact Maxwell's Fishpond. 
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FIG. 4: Snapshots from a Maxwell's Fishpond 200mm di- 
ameter and 10mm deep, simulated using OpenFOAM, at the 
start and each subsequent reformation up to the fifth. Al- 
though the simulation retains to a reasonable extent the re- 
peated refocussing, the wave dispersion continually increases 
the length of the pulse of ripples, so that reformations, while 
still remaining on the scale of a wavelength, have an ever in- 
creasing duration. 



Unfortunately, it is hard to build a real water wave 
device that works perfectly in the shallow water limit. 
Most notably, we will expect to see some residual dis- 
persion, due to the depth-dependent speeds of different 
wavelength ripples: the shallow water wave speed has 
some wavelength dependent speed variation (dispersion), 
it only matches eqn. @ for large wavelengths A. A more 
accurate wave speed equation is 



2 # A + u 
v = — tann 

2tt 



2nd 



(12) 



but even this neglects the effect of surface tension, which 
is discussed later in section Thus the brief reforma- 
tions expected in the ideal case will blur out into longer 
process as different wave components refocus at different 
times, and for the shorter A waves, the depth profile will 
be less perfect. 

As an estimator of the necesarily imperfect fishpond 
experiment, we (PK) also did more realistic finite ele- 
ment simulations using the open source simulator Open- 
FOAM (2(J using the inter MixingFoam engine, on a fast 
desktop PC. These simulations used a customised block 
meshing for the fishpond, because we wished to specify 
the depth profile at high resolution. Due to computa- 
tional constraints, the mesh could not be made as fine 
as we would like, but the general character of the ide- 
alised process was preserved, and the effect of dispersion 
demonstrated. A typical simulation result is shown on 
fig. ^ others indicate that larger fishponds may perform 
better than ours - although will be harder to construct, 
and will suffer more from dispersion. 




Screen 



FIG. 5: The experimental setup of NK and TT. The water 
dropper was held in place by a retort stand (not shown). 



IV. EXPERIMENTS 

To make the Maxwell's Fishpond, our departmental 
Mechanical Workshop machined us a brass insert with 
the necessary depth profile (see fig. [|) and mounted this 
in a nylon ring, with a small notch to indicate the pre- 
ferred water level. The choice of brass and nylon was 
based on convenience, not necessity, any waterproof ma- 
terials could suffice. Our brass insert was not perfectly 
smooth, but was machined to tolerances of well less than 
1mm. 

Two schemes for obtaining quantitative data were 
used: one by NK and TT, the other by CT and JT. Both 
were influenced by the fact that although viewing the 
fishpond directly gives a very strong impression of how 
well it works, this human perception does not translate 
easily into objective experimental data. Although even 
tiny ripples were surprisingly visible to the eye when in 
motion, it was less easy to get good experimental images. 

Being the first attempt, the NK/TT experimental 
setup was relatively simple. In a darkroom, they reflected 
lamplight off the water surface onto a screen, and took 
video images of the screen (see fig. |J). This screen was 
shielded from any direct light from the source lamp. The 
rippled water surface caused intensity variations on the 
screen, depending on whether the particular perturbation 
tended to focus or defocus the light (see fig. ^). Looking 
at light reflected from the ripples mean that the light ray 
paths diverged more strongly and made even very small 
ripples visible. The intensity pattern then indicated the 
progress of the ripples from source to image, and back. 
They then analysed the video by eye, frame by frame, to 
locate the times and positions of the reformations. 

The CT/JT setup did not involve a dark room, as one 
was not available at the time, instead they used lightproof 
box to eliminate stray light. Inside the box they imaged 
the reflection of a diffuse light source directly using a 
high resolution webcam (see fig. ^). This enabled them 
to electronically process the images. 

In both experimental setups, ripples were created using 
a water dropper to drop a single water droplet into the 
fishpond; the height of which was varied to ensure that 
the strongest possible ripples were generated, but not so 
big as to be accompanied by splashing or bubbles. A 
full range of starting positions was investigated, since 
the Fishpond should create reformations from any point 
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FIG. 6: A snapshot of the NK/TT experiment in progress, 
with the fishpond in the foreground, and the screen above and 
behind. We can clearly see the reflected ripple patten on the 
right of the screen; this image was taken just before its second 
reconvergence. 



Diffuse 
light 
source 




¥ Water 
dropper 



Webcam 



Light-proof box 

FIG. 7: The experimental setup of CT and JT. In addition 
the Fishpond was placed in a waterbath with a heat pump, to 
enable the temperature to be changed in a controllable way. 



- although starting positions near the wall suffered due 
to edge effects. 

The positioning of screens, light sources, cameras, and 
so on were systematically varied to achieve the best im- 
ages. For NK/TT, a shallow angle of reflection enhanced 
the images, although if too shallow this significantly re- 
duced the fraction of the water surface that could be 
seen. For CT/JT, the diffuse light source was placed at 
a low angle, but with the camera at a 90° reflectance an- 
gle; thus giving a strong contrast from the light reflected 
off the ripples. The diffuse light source avoided prob- 
lems caused by reflections off the bottom centre of the 
fishpond. 




FIG. 8: Typical differential image take from the webcam 
data of CT/JT, after speckle was removed. The purple lines 
cross the centre of the fishpond, with one linking source-point 
to reformation point, and the other at right angles. The image 
is foreshortened due to the angle of view of the webcam. 



A. Setup 1 

Using the setup of fig. ||, videos were taken from each 
initial droplet to final dissipation of the ripples. Each 
instance was viewed carefully frame-by-frame to deter- 
mine the time taken for the 1st and 2nd reformations. It 
should be noted that reformations occuring after the sec- 
ond were not quantitatively analysed because their rip- 
ples were too faint to measure with an acceptable error 
in accuracy - although they were still visible to the eye 
as a fleeting event. 

Reformation times were found by first locating the de- 
sired reformation within a short "target" sequence of 
video frames. With each reformation point taken to be 
indicated by the frame with the most localised ripple pat- 
tern, the short target sequence was viewed independently 
by NK and TT, both forward and backwards. Once the 
best reformation frame (s) were chosen, the frame indexes 
were converted into a time using the frame rate of the 
camera, and averaged if necessary. The times from the 
start to the first and second reformations could then be 
calculated. 

The positions of source, and 1st and 2nd reformation 
points were also compared and measured. However, since 
the foreshortening of the reflected ripple pattern means 
that it appears as an ellipse, times were calculated from 
instances where the source and reformation points ap- 
peared on a horizontal line on the screen, i.e. the long 
axis of the ellipse. 

The average times taken for 1st and 2nd reformations 
were (1.04±0.04)s and (1.07 ± 0.12)s. The variation be- 
tween times taken for 1st reformations of different ripples 
at varied start positions was found to be 0.04s; and so 
is dominated by the discretization of the time intervals 
between each frame of video. The variation for the sec- 
ond reformations is larger, and was largely due to the 
difficulty in determining the video frame with the most 
localized ripple pattern, since by then the ripples had not 
only suffered increased dispersion but also greater dissi- 
pation. 

The theoretically expected reformation time can be 
calculated most easily by referring to the reference 
sphere. The water depth at the centre of the Fishpond 
matches that on the imaginary reference sphere, which 
in our case is 2.5mm deep, leading to a wave velocity of 
v = \J~gd ~ 157mm/s. Each reformation on the sphere 
take place after one half circumnavigation, and the sphere 
radius s is half that of the Fishpond radius of ro = 10cm; 
thus the theoretical reformation time is 



_ tts 3.142 . 50mm 



157mm/s 



1.003s, 



(13) 



which seems to be in good agreement with the measure- 
ments - but see also the discussion in section 

NK and TT also attempted to increase the longevity 
of the ripples, to see if more reformations could be mea- 
sured. In one trial, they tried reducing the water viscos- 
ity (and hence dissipation) by using water at 70 C in the 
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Fishpond, but this also had negligible effect. In another, 
they tried to affect the surface tension using both deter- 
gent and temperature (separately), but this slowed the 
wave speeds, so that more dissipation occurred between 
each reformation. 



B. Setup 2 

Rather than analysing video of ripple reflections frame 
by frame, CT and JT instead directly imaged the ripples 
using a webcam and VirtualDub |2l| , using the setup of 
fig. [7[ It was helpful that as the video was being im- 
aged, the software also reported on the framerate, which 
helps in converting frame indexes into times. Next, to 
bring out the ripple dynamics from the background im- 
age, they subtracted each frame from the previous one 
using AviSynth [Q. The resulting differential image 
data as shown in fig. ||, reveals only the wave motion, 
which was then analysed using Tracker |^3| . Each image 
was scanned along the axis between source and image 
points, giving a 2D dataset, comprising a time series of 
ID datasets along this axis. The result was processed an 
plotted using a variety of software, including Microsoft 
Excel, Matlab §|, and Scilab f|. 

In order to optimise the reformation process, CT and 
JT systematically analysed results taken for a range of 
temperatures and fill volumes, as shown in tables | and 
|J. The optimum temperature was found to be 15C: al- 
though the water viscocity (and loss) increases for lower 
temperatures, the surface tension increases, leading to 
faster wave speeds. Note that the optimum fill volume 
centred around 190ml, this can be compared to that for 
the design parameters, which radial integration of the de- 
sign depth profile is found to be 183ml. This is in good 
agreement - a 3ml change in fill volume corresponds to 
about a 0.1mm depth change, and velocity shifts of be- 
tween 2% or less. The experimental optimum filling vol- 
ume of 190ml is higher than the design volume, perhaps 
because of the way the Fishpond fills - e.g. the design 
takes no account of surface tension. 

In fig. [9] we see the first reformation very clearly, al- 
though dispersion has spread out the initial impulse in- 
troduce by a falling water drop. On the first traversal, 
we can not only see the evidence of several ripple crests, 
but also the slight fanning as the ripples disperse. Fur- 
ther reformations, although apparent to the eye, do not 
show up over the imaging noise. 



Volume [ml] 


170 


175 


180 


185 


190 


195 


200 


Reformations 


2 


2 


2 


3 


3 


3 


2 



TABLE I: Relationship between water volume in the fishpond 
and the number of reformations visible to the eye. 

Reformation times and errors were extracted from data 
like that shown on fig. fusing a curve fitting process. For 
example, at the optimum temperature of 15C, the first 




1 1 1 1 

10 30 50 70 

Frame Index 

FIG. 9: Differential luminance data as a function of time 
taken along the axis between initial disturbance and refocus 
point, for two different water temperatures: (a) 15C, and (b) 
20C. These contour plots are made by averaging over adjacent 
points and using a logarithmic scale; the contours are evenly 
spaced, with a minimum level chosen to best display the ripple 
patterns. 



Temperature [C] 





5 


10 


15 


20 


25 


30 


35 


40 


Reformations 


1 


2 


2 


3 


2 


2 


2 


2 


2 



TABLE II: Relationship between temperature of the fishpond 
and the number of reformations visible to the eye. 



recurrence was calculated to have occured at the 23rd 
frame (at 0.72 ± 0.06s), with the second being 30 frames 
later (+ 1.03 ± 0.11s). 

As noted already, the recurrence time including sur- 
face tension effects is faster than that of about Is calcu- 
lated in eqn. ( |l3| ) by some 30%; in quite good agreement 
with the measured value. The second recurrence was 
notably slower than the first, but then the error is also 
much larger; but of course some slowing might be ex- 
pected since the longer wavelengths both persist longer 
and travel more slowly. 

Finally, despite the difficulty in extracting second and 
third reformation times from the video data, and in see- 
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ing them in plots such as fig. |9|, in the original videos 
themselves the third reformation is clearly visible to the 
eye. This suggests that significant performance improve- 
ments are still possible in the automated processing of 
the data. 



V. SURFACE TENSION 

Since water is the obvious liquid to use in the Fish- 
pond, and it has a significant surface tension, we should 
estimate its effects. 

The wave velocity including the effects of surface ten- 
sion (7 (in N/m) on waves in a fluid with densi ty p and 
wavelength A = 2n/k is not a simple as in eqn. (II4), but 
is 



k? 



1 



gp 



tanh (kh) . 



(14) 



Thus the correction to the leading term which gives us the 
shallow water wave speed is a factor e = ak 2 /gp. Surface 
tension in water reduces with temperature, being about 
0.073 N/m at 20C, but 0.061 N/m at 90C. At about 20C, 



ak 2 ^ 0.073/c 2 
~gp ~ 9.81 x 1000 



= 7.44 x 10" 6 /c 2 . 



(15) 



For water waves of wavelength 2cm, k = 27r/0.02 « 
314m _1 , so that e w 0.73; the wave speed is therefore 
a factor of \/1.73 (or 30%) higher than expected based 
on depth alone, and is wavelength dependent even in the 
shallow water limit. However, this speed shift is not depth 
dependent, so the refocussing character of the Fishpond 
is unaffacted - but different wavelengths reform at differ- 
ent times. This dispersion means that a determination of 
the reformation time becomes harder as time progesses, 
and will depend on the specific details of how a given 
reformation time is evaluated. 

Initially, the NK/TT measured reformation time of 
(1.04 ± 0.04)s seemed i n g ood agreement with the sim- 
ple prediction of eqn. (|l3|), i.e. 1.00s. However, we can 
now see that surface tension effects should reduce this 
prediction by 30% to about 0.7s, which agrees with that 
measured by CT/JT, and not that of NK/TT. 

But why do the two experiments give such different 
outcomes? Two scenarios, which are not mutually exclu- 
sive, suggest themselves. 

First, the criteria for choosing the reformation times 
differed, and this will affect which frame of video selected 
- NK/TT chose by eye the frame with the smallest re- 
gion of disturbed water, whereas CT/JT applied a simple 
fitting algorithm to digitised data along one axis. 

Second, NK/TT relied upon the reported frame rate 
of their camera, and perhaps this was not reliable; al- 
though with hindsight we realise that their framerates 



might have been easily calibrated by videoing a clock ei- 
ther before or at the same time as the each experimental 
run. 

It is gratifying that the more sophisticated setup of 
CT/JT gives good agreement with theory, although it is 
not clear why the first attempt by NK/TT did less well. 
Neverthless, one of the features of this student project 
was it could be implemented in many different ways - the 
students were given the Fishpond, some reading material, 
and some suggestions and then largely left to get on with 
it independently. Still other experimental set-ups and 
measurements are possible, and so we expect that the 
Fishpond itself will be reused many times in the future. 



VI. SUMMARY 

We have shown how an exotic phenomenon from trans- 
formation optics - the Maxwell's Fisheye lens - can 
be converted into simple water waves in a tabletop 
"Maxwell's Fishpond". This is currently being used 
sucessfully as a third year undergraduate experimental 
project in the Physics Department of Imperial College 
London. While the remarkable series of image reforma- 
tions provides the hook which makes the project inter- 
esting, there are many other features that can be investi- 
gated as part of the experiment. Most straightforwardly, 
there is a variety of imaging possibilities to be investi- 
gated (two of which were discussed here), and various 
experimental conditions - lighting, fill depth, etc - to 
be determined. Also, the effect of viscosity on perfor- 
mance can be tested by changing the water temperature, 
or surface tension can be removed by adding detergent - 
or other liquids might be used. A transparent Fishpond 
might be made so as to image the ripples in transmis- 
sion, or a vibrating source could be use in an attempt to 
generate standing waves. For the more mathematically 
inclined, the nature of stereographic projections can be 
researched, other comparable devices - e.g. the Eaton 
or Luneburg lenses - considered, or numerical simula- 
tions attempted. Alternatively, rather than only aiming 
to optimise the number of reformations, or visibility to 
the eye, but it is also possible to consider ease or sim- 
plicity of fabrication, with a view to testing performance 
as a function of size. Since our simulations show that 
the general behaviour persists even for an approximate 
depth profile - such as a shallow dome - sophisticated or 
precise manufacturing processes are not needed. 

In this way, this simple, eye-catching device provides a 
rich playground in which a wide variety of students can 
test their skills while investigating a novel device not only 
part of contemporary research - that of transformation 
optics and acoustics - but with a history that goes back 
to Maxwell himself. 
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